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Sood and Grassberger studied in |lj random walks on 
random graphs that are biased towards a fixed target 
point. They put forward a critical bias strength be such 
that a random walker on an infinite graph eventually 
reaches the target with probability 1 when b > be, while 
a finite fraction of walks drift off to infinity for b < be- 
They rely on rigorous results obtained for biased walks on 
Galton- Watson (GW) trees [1 to calculate be, and give 
arguments indicating that this result should also hold 
for random graphs such as Erdos-Renyi (ER) graphs and 
MoUoy-Reed (MR) graphs. To validate their prediction, 
they show by numerical simulations that the mean return 
time (MRT) on a finite ER graph, as a function of the 
graph size N, exhibits a transition around the expected 
be- 

Here we show that the MRT on a GW tree can actually 
be computed analytically. This allows us (i) to show 
analytically that indeed the MRT displays a transition 
at be, (ii) to elucidate the N dependence of the MRT, 
which contradicts the oc N scaling expected in [l| for 
b < be. 

Let us consider a realization of a GW tree of g genera- 
tions, and denote by z„ the number of nodes of generation 
n (0 < n < g — 1). The number of nodes at generation 
zero is taken equal to one [zq = 1) and each node has 
a random number of daughter nodes of mean k. As in 
ref. we consider a random walker starting from the 
root node {n = 0) and experiencing a constant bias, such 
that the probabilities pf and to jump from site I re- 
spectively towards and away from the root are given by 
pj' = b/Afi and = b~^ /Mi where A/; is a normalization 
constant. 

The key point of the derivation is that the MRT, here 
denoted by (T) , is given by the Kac formula 0, 0, Bl : 
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tion n which is easily showed to verify Pcq(n) cx Znb 
Normalization then yields straightforwardly 
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We then denote by X the average of a quantity X over 
the realizations of the graph. Using = fc" we 
first obtain the average number of nodes in the GW tree 
= (1 — — k), and finally the desired quantity 
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with e = ln(fc/62)/ln(fc). 

Equation ([3]) clearly shows that (T) exhibits a transi- 
tion when e = or equivalently b — be — \fk, in agree- 
ment with the numerical simulations of 1]. Furthermore, 
we obtain explicitly the asymptotics of (T) for large A^: 


(T) 


fc/b2-^l ^' for e > (6 < be) 
for e = (6 = 6,) 
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where Pca^(n) is the equilibrium distribution at genera- 


Equations ([3]) and ([4]) show that in the case of an unbi- 
ased walk (6 = 1), one has (T) — N, but that this scaling 
does not hold as soon as 5 7^ 1. Note however that this 
exact scaling could be hard to distinghuish from the nu- 
merical simulations of 

We thank the authors of [l| and acknowledge that they 
have agreed on the validity of this result and agreed to 
publication of this Comment. 
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